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ABSTRACT 

Recent infrared (IR) observations of freshly-formed dust in supernova remnants (SNRs) have yielded 
significantly lower dust masses than predicted by theoretical models and measured from high redshift 
observations. The Crab Nebula's pulsar wind is thought to be sweeping up freshly-formed supernova 
(SN) dust along with the ejected gas. The evidence for this dust was found in the form of an IR excess 
in the integrated spectrum of the Crab and in extinction against the synchrotron nebula that revealed 
the presence of dust in the filament cores. We present the first spatially resolved emission spectra of 
dust in the Crab Nebula acquired with the Infrared Spectrograph aboard the Spitzer Space Telescope. 
The IR spectra are dominated by synchrotron emission and show forbidden line emission from from S, 
Si, Ne, Ar, O, Fe, and Ni. We derived a synchrotron spectral map from the 3.6 and 4.5 /xm images, and 
subtracted this contribution from our data to produce a map of the residual continuum emission from 
dust. The dust emission appears to be concentrated along the ejecta filaments and is well described 
by an amorphous carbon or silicate grain compositions. We find a dust temperature of 55 ± 4 K for 
silicates and 60 ± 7 K for carbon grains. The total estimated dust mass is (1.2 — 12) x 10~ 3 M , 
well below the theoretical dust yield predicted for a core-collapse supernova. Our grain heating model 
implies that the dust grain radii are relatively small, unlike what is expected for dust grains formed 
in a Type IIP SN. 

Subject headings: dust, extinction - infrared: ISM - ISM: individual objects (Crab Nebula) - ISM: 
supernova remnants - pulsars: individual (PSR B0531+21) 



1. INTRODUCTION 

Supernova (SN) dust is expected to form and has been 
observed to form shortly after the SN explosion in regions 
where the cooling gas is dense enough fo r grain growth 
to occur (see review bv lGall et afll20Tl . The spectral 
signatures of this SN-formed dust depend on its temper- 
ature and composition, properties that are determined 
by the available gas abundances and the environmental 
conditions in the supernova remnant (SNR). Theoretical 
models based on classical nucleation theory predict that 
a range of dust species can form in SN ejecta, with some 
of the most abundant being Mg SiOs, SiC>2, Mg 2 Si04, Si, 
and C (e.g.. iKozasa et aLll2009f ). The same models pre- 
dict that the total dust mass produced per SN explosion 
should be 0.1-0.7 M Q , with 0.04-0.2 M Q surviving the 
reverse shock and being recycled back into the interstel- 
lar m edium (ISM) fe.g.. iDwek et al.ll2008t IKozasa et all 
2009). The large quantities of dust observed in high- 
redshift galaxies could be e xplained if an av erage SN 
produced 0.1-1.0 M (e.g.. IDwek et alj|2009f) . but re- 
cent IR observations of dust in SNRs have yielded much 
lower masses. 
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Detailed studies of dust emission spectra from SN- 
formed dust have been limited to the closest and 
youngest SNRs, mainly due to the lack of spatial res- 
olution that is required to spatially separate ejecta dust 
from the swept-up ISM dust and confirm its SN origin. 
In recent years, spectroscopic observations of supernova 
remnants (SNRs) with the Spitzer Space Telescope have 
confirmed the ejecta origin of SN dust in a few rem- 
nant s by correlating its spatial distribution with ejecta 
gas (|Rho e t al. 2009]) . The observed emission spectra re- 
vealed spectral features that differ from those of typical 
ISM dust that predominantly consists of astronomical sil- 
icates and carbonaceous material, such as polycyclic aro- 
matic hydrocarbons (PA Hs), graphite, and amorphous 
carbon (for a review see iZubko et all 120041 ) . IRho et al.l 
(2009) identified ejecta dust in the SNRs 1E102.2-7219, 
N132D, and Cas A, and found emission features from 
several dust species, including Si02, MgSiOa, Si, SiC, 
AI2O3, Fe, and featureless carbon dust. IR emission 
spectra of ejecta dust in Cas A and SNR G54.1+0.3 
revealed a b road 21 /xm dus t feature that still remains 
unidentified (|Rho et al.ll200l iTemim et al.l[2010h . 

The question of whether SNe are major contribu- 
tors of dust to the ISM has still not been settled. 
The mass estimates of newly formed warm dust from 
recent Spitzer observat i ons a r e in the 0.02-0 . 1 Mq 
range (ISugerman et al.l [2001 IRho et al.l [200l [2001 
Temim et all 120101: IBarlow et al.l 120101: ISibthorpe et al.l 
20101). Herschel obs ervations of SN 1987A revealed 0.4- 
0.7 M of cool dust ([Matsuura et al.ll201lD . while Spitzer 
observations reveale d up to 0.1 M^ of ne wly formed dust 
in SNR G54.1+0.3 (jTemim et alJl20Tot) . However, this 
dust will eventually encounter the SNR reverse shock. 
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Without detailed information about the grain composi- 
tion and mass, it is difficult to estimate how much of this 
dust will survive the encounter and be injected back to 
the ISM. 

The Crab Nebula is a prototypical example of a pul- 
sar wind nebula (PWN) that is sweeping up SN ejecta 
material, observed in the form of brig ht optical filaments 
(jGaensler fc Sland [20061: lHester|[200l . The evidence for 
the existence of dust in the Crab Nebula was found in 
the form of an IR ex cess above the PWN's synchrotron 
power- law spectrum (fR imblc 19 771: iGlaccum et al.lll982t 
Marsden etahl 119841: iDouvion et al.l 120011 : iGreen et al.l 



20041: iTemim et al.ll2006ft . The IR excess observed in the 



integrated spectral energy distribution (SED) of the Crab 
implies a total dust mass of 0.001-0.03 M Q , a range that 
depends on assumptions about dust composition and 
temperature. Optical observations revealed absorption 
features from dust that sp atially correlate with the cores 
of the optical filaments dWoltier fc Veron-Cettvl 
Fesen fc Blahlll990l; IBester et al.lll990l: iBlair et all 



1987 



1997 



Loilll2010[) . Based on the measured dust-to-gas ratios in 
the filaments that range from 0.03-0.1 and the spatial 
correlation between dust extinction and low ionization 
lines, the dust in the Crab is most li kely SN-formed dust 
that has b een swept up by the PWN (|Sankrit et al.lll998l: 
lloll[20rol) . 

The properties of the dust in the Crab Nebula are not 
well determined. The IR continuum is dominated by syn- 
chrotron emission, and prior to Spitzer, there were no 
detailed detections of dust emission spectra that could 
shed light on the grain composition. In this work, we 
removed the synchrotron contribution from spatially re- 
solved Spitzer spectra of the Crab Nebula in order to 
isolate the dust emission spectrum for the first time and 
determine the physical properties and spatial distribu- 
tion of the SN dust. 

2. OBSERVATIONS AND DATA REDUCTION 

The IR spectroscopy was carried out at six differ- 
ent positions across the Crab Nebula with all four 
modu l es of the Spitzer Space Telescope's (jWerner et al.1 
2001 IGehrz et al l [200l Infrared Spectrograph (IRS, 
Houck et all 120041) : short-high (SH, 9.9-19.6 ,um), long- 
high (LH, 18.-9-37.2 jum), short-low (SL, 5.2-14.3 /mi), 
and long-low (LL, 14.1-38.0 /im). In this paper, we 
present sample SL and LL spectra. The observations 
were carried out between 2004 and 2006 under project ID 
24, and processed with the pipeline version S18.7.0. They 
are summarized in Table [TJ including the coordinates, ex- 
posure times and AORKEYs. The corresponding slit po- 
sitions are shown in Figure [TJ A background spectrum 
was obtained with all four IRS modules at a position 
approximately 5. '4 southwest of the center of the nebula 
(see Table[TJ). The spectra were extracted, calibrated and 
cleaned using the CUbe Builder for IRS Spectral Maps 
(CUBISM) version vl.7 ([Smith et al.l2007l ) and the back- 
ground spectra for each module were subtracted from 
the source spectra. The background does not appear to 
vary spatially in the vicinity of the Crab Nebula, nor 
does it contain any strong emission lines that would con- 
taminate the spectrum. Depending on the wavelength, 
the continuum level in the source spectrum is 3-7 times 
higher than the background level. The analysis was car- 
ried out with Spectroscopic Modeling, Analysis and Re- 



FlG. 1. — Positions of the IRS slits overlaid on the MIPS 24 fim 
image. The LL slits are the large slits oriented in the East/West 
direction and span the entire length of the SNR shell, and the SL 
slits are the narrower, North/South oriented slits. 



duction Tool (SMART) version 8.1.2 (Hig don et al.|[20p: 
iLebouteiller et al.ll201ll) . 

The Spitzer Infrared Array Camera (IRAC) observa- 
tions at 3.6, 4.5, 5.8, and 8.0 /zm were carried out on 
2004 March 6, while the Multiband Imagining Photome- 
ter for Spitzer (MIPS) observations at 24 and 70 /im were 
carried out on 2004 March 14. The details about the re- 
duction of th e IRA C and MIPS data can be found in 
ITemim et al.l (|2006[ ). where the images were previously 
published. All images in this paper are oriented such 
that North is up and East to the left. 

3. ANALYSIS AND RESULTS 

3.1. Subtraction of the Synchrotron Contribution 

The synchrotron emission from the Crab Nebula can 
be described by a power law spectrum of the form 
L y oc v~ a , where L v is the synchrotron specific luminos- 
ity and a is the spectral index. At radio wavelengths, 
the s pectral index is given by a — 0.3 (Baa rs et al.l 
I1977D. and steepens to a = .7 at optical wavelengths 
(jVeron-Cettv fc Woltierl I1993D . This suggests that a 
break in the power-law spectrum occurs somewhere i n 
the IR, between 10 and 1000 /im (|Marsden et al.lll984(K 
The most recent estimate for this break frequency is 



(7.0 



+6.5\ 
-3AJ 



10 2 GHz, or (4.3+2:1) x 10 A* m > wel1 beyond 



the Spitzer wavelength range (|Arendt et al.H201lD . 

The first step in estimating the synchrotron contribu- 
tion at each spatial position in the Crab Nebula was to 
create a map of the synchrotron spectral index a, We 
used the synchrotron-dominated IRAC 3.6 /im and 4.5 
/im images to compute an index map, given by the equa- 
tion 

ln(F 1 /F 2 ) 



a 



/n(Ai/A 2 ) 



where F is the surface brightness and A is the wave- 
length. The IRAC images were first background sub- 
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TABLE 1 
IRS OBSERVATIONS 



Position 


RA/Doc (J2000) 


AORKEY 


Module 


Exposure/Cycle (s) 


Cycles 


Observation Date 


Position 1 


5i 1 34™31.00?+22°01'10'.'0 


3857664 


SL (LL) 


480 (240) 


2 (4) 


2004-03-05 


Position 2 


5! 1 347 1 29.5ffi+22 o 00'31<'0 


16200704 


SL (LL) 


60 (30) 


8 (8) 


2006-03-07 


Position 3 


5! 1 347 1 33.88?+22°00'52<'0 


16201216 


SL (LL) 


60 (30) 


8 (8) 


2006-03-07 


Position 4 


34™ 34. 53? +2 1 ° 59' 46'.' 5 


12634624 


SL (LL) 


60 (30) 


8 (8) 


2005-10-13 


Position 5 


5^34™38.31?+22°00'49'.'2 


12634624 


SL (LL) 


60 (30) 


8 (8) 


2005-10-13 


Position 6 


5^34™38.44?+21°59'l0''0 


16200704 


SL (LL) 


60 (30) 


8 (8) 


2006-03-07 


Background 


5! 1 34™ 39.91! +22°05'59?6 


16200704 


SL (LL) 


60 (30) 


8 (8) 


2006-03-07 




5^34™39.91?+22°05'59('6 


12634624 


SL (LL) 


60 (30) 


8 (8) 


2005-10-13 



Note. — The spectra at position 1 were previously published by Tcmim ct al. (2006). 



tracted and corrected for extinc t ion b ased on the A\j Ak 
values from llndebetouw et al.l (|2005l), and the re lation 
N H /A K = 1.821 x 10 22 cm" 2 from |Draine] (fl989t) . The 
hydrogen column density for the Crab wa s set to Nh = 
3.54 x 10 21 cm" 2 (jWillingale et al.ll200l . This led to 
correction factors of 1.11 and 1.08 for the 3.6 and 4.5 /im 
channel, respectively. The images were also corrected 
by the IRAC surface brightness correction factors for ex- 
tended sources; 0.91 for 3.6 /im and 0.94 for 4.5 /im. The 
3.6 /im image was then convolved to the resolution of the 
4.5 /im image using t he convolution kernels developed by 
iGordon et al.l (|2008| ). The nominal wavelengths for the 
two channels that were used in the final calculation for 
the index are 3.550 /im and 4.493 /im, as specified by the 
IRAC Instrument Handbook. 

The final synchrotron index map is shown in Fig- 
ure [21 Due to the updated extinction correction, the 
va lues of the index ar e slightly higher than those found 
in iTemim et al.l (|2006T ) . The value of the synchrotron in- 
dex in the inner torus and jet is approximately a = 0.3 
and it gradually steepens towards the outer edges of the 
nebula due to synchrotron losses. The index steepens to 
about a = 0.6 towards the SE, and a = 1.0 towards the 
NW edge of the nebula. A similar asymmetry is seen 
in the spectral index map at optical wavelengths, where 
the index steepens to a ~ 0.8 in the SE, and e xceeds 
a = 1.0 in the NW (IVeron-Cettv fc Woltierl[l993l) . The 
calibration uncertainty on the ratio of the fluxes in the 
two I RAC channels is approximately 4% (iReach et al.l 
2005), which translates into a ±0.18 uncertainty in the 
synchrotron spectral index. This large uncertainty in the 
spectral index is systematic across the nebula and does 
not affect the spatial variations in the in index. 

In order to investigate the morphology of the dust and 
line emission in the Crab Nebula, we subtracted the syn- 
chrotron contribution from the IRAC 5.8 and 8.0 /im im- 
ages and the MIPS 24 and 70 /im images. We produced 
a synchrotron-subtracted image for each wavelength by 
subtracting the synchrotron contribution extrapolated 
from the IRAC 3.6 /im image using the equation 

Fi,residual = Fi ~ ^3.6(3.550/A 2 ) Q . (2) 

Here, the subscript "i" refers to the wavelength of the 
image from which the extrapolated synchrotron emission 
is being subtracted. In each case, t he 3.6 /im image wa s 
convolved to the resolution of A^ (jGordon et al.l 12008). 
and each of the ima ges was extinction correc ted using the 
extinction curve of IChiar fc Tielensl (|2006t) . The IRAC 



surface brightness correction factors were applied to all 
IRAC images (Table 4.9 of the IRAC Instrument Hand- 
book) . The resulting synchrotron-subtracted images are 
shown in Figure [3] and will be discussed in Section |H 

Next, we subtracted the synchrotron contribution from 
the SL and LL spectra for each spatial position. We cre- 
ated four synchrotron spectral cubes covering the wave- 
length range of each of the four low resolution IRS mod- 
ules for every spatial pixel along the slits shown in Fig- 
ured We first checked the calibration of the IRS spectra 
by comparing the LL emission integrated over the MIPS 
24 /xm bandpass to the flux measured in the MIPS 24 
/im image. We found the normalization factor between 
the LL spectra and the 24 /im surface brightness to be 
0.997. There was no need for an additional scaling fac- 
tor for the SL spectra, since the continuum levels in the 
overlapping wavelength regions of the LL and SL spectra 
matched well. We also integrated this final spectrum over 
the IRAC 8.0 /im bandpass and found that the flux does 
not differ significantly (< 5%) from the flux in the IRAC 
8 /im image. We mapped the 3.6 /im image to match 
the spectral cube map for each of SL and LL orders, and 
used equation [5] to compute a synchrotron emission cube 
for each module. These synchrotron cubes were then 
subtracted from the background-subtracted IRS spectral 
cubes that have been extinction correcte d using the ex- 
tinction curve of IChiar fc Tielensl ((2006). This resulted 
in SL and LL synchrotron-subtracted spectral cubes that 
only include line emission and continuum emission from 
dust. 

3.2. Residual Dust Emission 

Inspection of the synchrotron-subtracted spectral cube 
reveals residual continuum emission from dust in some 
regions of the nebula. In order to determine the spa- 
tial distribution of the dust emission in the IRS slits, we 
used the synchrotron-subtracted cube to produce a map 
of the rising continuum. We selected a line-free region 
of the spectrum and produced a map of the integrated 
emission between 27-32 /im. This map is shown in Fig- 
ure H] in magenta, overlaid on the MIPS 24 /im image in 
blue. The map shows that the brightest dust emission co- 
incides with the brightest filaments in the 24 /im image. 
We note that the background continuum in the vicinity 
of the Crab Nebula falls off longward of 25 /im, so the 
background contamination in the excess dust continuum 
is minimal. 

To extract spectra for analysis of the dust emission, 
we chose three positions on the dust map that provided 
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TABLE 2 
IRS LINE FITS 



Line ID 


— 7T, — ; ; — 

Wavelength 
(H 


Region 1 Flux 
(lO" 11 W/cm 2 /sr) 


Hcgion 2 b lux 

f -\ r\ 1 1 tit / 9/ \ 

(10 11 W/ cmr/sr) 


Hegion 6 b lux 
(10 u 'W/cm z /sr) 




Fe III 


5.340 


3.897 


± 


0.299 


1 937 + 239 


1 579 + 174 




Me Vl 


5.610 








n 929 + 040 






Ni II] 


6.636 


2.474 


± 


0.170 


948 4- 099 


712 + 062 




Ar II] 


6.985 


16.72 


± 


0.457 


q 7KQ _|_ o i sn 


O . VJ O C< _1_ VJ.-LiJV^ 




Na III] 


7.318 


0.365 


± 


0.086 








Ar III] 


8.991 


3.308 


± 


0.007 


2 093 + 526 


1 481 + 039 




s iv] 


10.511 


2.372 


± 


0.204 


2 256 + 100 

£i,£nJ\J _1_ U.1UU 


i ™q 4- n nc;t; 

-L.iJC'iJ _1_ Kl • \J <J <J 




Ni II] 


10.682 


0.222 


± 


0.004 


n ni o + o oo6 






Ni II] 


11.308 


0.067 


± 


0.007 








Ne II] 


12.814 


9.383 


± 


0.350 


46.05 ± 1.14 


19.32 ± 0.41 




Ne V]/[C1 II] 


14.322/14.368 


0.956 


± 


0.220 


4.320 ± 1.501 


1.831 ± 0.299 




Ne III] 


15.555 


7.565 


± 


0.123 


34.32 ± 0.038 


18.77 ± 0.318 




Fe II] 


17.936 


0.560 


± 


0.010 


0.491 ± 0.011 


0.441 ± 0.027 




S III] 


18.713 


4.278 


± 


0.042 


4.064 ± 0.137 


4.895 ± 0.308 




Fe III] 


22.925 


0.322 


± 


0.024 


0.524 ± 0.046 


0.268 ± 0.037 




Ne V]/[Fe II] 


24.318/24.519 


0.281 


± 


0.034 


1.215 ± 0.087 


0.924 ± 0.047 




O IV] /[Fe II] 


25.890/25.988 


3.913 


± 


0.010 


11.08 ± 0.34 


8.739 ± 0.278 




S III] 


33.481 


3.876 


± 


0.112 


3.321 ± 0.097 


4.194 ± 0.126 




Si II] 


34.815 


2.684 


± 


0.151 


2.400 ± 0.125 


3.057 ± 0.143 




Fe II] 


35.349 


0.243 


± 


0.035 








Ne III] 


36.014 


0.510 


± 


0.014 


2.397 ± 0.134 


1.279 ± 0.117 



Note. — The listed uncertainties are 1-cr statistical uncertainties from the fit only. 




0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

Fig. 2. — Synchrotron index map calculated from the 
synchrotron-dominated IRAC 3.6 and 4.5 (im images. The spec- 
tral index a steepens from 0.3 in the torus and jet structures to 
0.7 (1.0) in the outer SE (NW) parts of the nebula. 

the best signal-to-noise. These regions are shown in Fig- 
ure 21 and also correspond to positions 2, 4, and 5 in 
Table Q] While the LL spectra were extracted from re- 
gions outlined in black, the SL spectra were extracted 
from smaller sub-regions that overlap with the LL slits 
(see Figure [IJ. In order to fit the dust emission only, we 
extracted spectra from the synchrotron-subtracted spec- 
tral cubes and measured the MIPS 70 /im flux density 
for these same regions. We fit the data with a modified 
blackbody model with three different grain compositions. 
The best-fit dust models and the unsubtracted spectra 
are shown in Figure \5[ The details of the fitting are dis- 
cussed in Section [5] 



3.3. Line Emission 

Strong forbidden emission lines are present in the IRS 
spectrum at all six positions in the Crab Nebula. We 
measured line intensities in the sample spectra from three 
positions in Figures 0] and [S] The measured line inten- 
sities are listed in Table [21 The strongest line emission 
comes from sulfur, neon, argon, oxygen and iron. Other 
relatively bright lines include silicon and nickel. There 
are some obvious differences in the relative line intensities 
between Region 1, the brightest filament in the south, 
and Regions 2 and 3, the brightest filaments in E/W 
direction. While the southern filament shows stronger 
emission from [Ar II], [Ar III], [Na III], and [Ni II], 
the east and west filaments show stronger emission from 
[Ne II], [Ne III] and [O IV]. All three regions show sim- 
ilarly bright emission from [S III], [S IV], [Si II], and 
[Fe II] . Emission from lines with higher ionization poten- 
tials, such as [O IV] , [Ne V] , and [Mg V] , is strongest in 
Region 2. Comparisons of line intensity ratios for lines 
of similar ionization potentials reveal some differences 
across the three regions that may be indicative of vari- 
ations in the gas abundances. While [S III] is equally 
bright in all three regions, [Ne II] is much more pro- 
nounced in Region 2. The [Ni II] /[Si II] line ratio is more 
than a factor of two higher in Region 1, and [Na III] is 
not even detected in Region 2, where [Ne III] and [O IV] 
are most prominent. 

Since the velocity splitting of the line emission from the 
two sides of the expanding nebula is not resolved in SL 
and LL modules, we cannot make any definitive conclu- 
sions about the variations in gas properties based on the 
low-resolution data. A detailed analysis and modeling 
of the high-resolution line emission detected by Spitzer 
will be presented in a subsequent publication. In this 
paper, we use the measured line intensities to determine 
the contribution of line emission to the total integrated 
flux from the Crab Nebula in order to estimate the total 
mass of dust. 
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Fig. 3. — Synchrotron-subtracted images of the Crab Nebula. The IRAC 5.8 and 8.0 fim images are on a linear gray scale from 0.0-10 
MJy/sr, and the MIPS 24 and 70 (im images are on a linear scale from 0-200 MJy/sr. The dashed circle indicates the region where the 
extrapolated synchrotron emission is overestimated and where the residuals are negative. 



4. MORPHOLOGY OF THE DUST AND LINE EMISSION 

The synchrotron-subtracted images of the Crab Neb- 
ula are shown in Figure El w hile the unsubtracte d im- 
ages are shown in Figure 3 of iTemim etaLl (|2006ft . The 
synchrotron component that was extrapolated from the 
IRAC 3.6 fim image describes the synchrotron emission 
at the longer IRAC and MIPS wavelengths fairly well. 
The synchrotron emission from the torus and jet features 
is completely subtracted out by our method. The ex- 
trapolation overestimates the synchrotron contribution 
in the NW part of the nebula, in a region shown by the 
dashed circle in Figure [3] The residual IRAC images 
show slightly negative values in the NW, and this over- 
estimate increases with wavelength, with the most sig- 
nificant over-subtraction at 70 /im. The over-subtracted 
values in the NW are on the order of ~ 1 MJy/sr at 5.8 
lira., 3 MJy/sr at 8.0 /im, 5-10 MJy/sr at 24 /im, and 15- 
25 MJy/sr at 70 /xm. These values range from approxi- 
mately 5-20 % of the total extrapolated synchrotron flux 
in this region, and do not appear to be caused by the un- 
certainties involved in the subtraction. 

The uncertainty in the spectral index is dominated by 



the calibration uncertainties in the IRAC flux density 
and extinction correction at 3.6 and 4.5 /zm, so it is ex- 
pected to be a systematic uncertainty with no signifi- 
cant spatial variations. The over-subtraction from the 
extrapolated synchrotron spectrum is seen in one local- 
ized region of the nebula, and this suggests that it may 
be caused by an intrinsic difference in the synchrotron 
spectrum in this region. As the magnetic field evolves 
in the expending nebula, we might expect a curvature 
in the PWN's power-law spectrum. A slight curvature of 
the spectrum in the IR may explain why the synchrotron 
over-subtraction increases with wavelength, but we still 
do not understand why the curvature may be different or 
more pronounced in NW part of the nebula in particular. 

The residual filamentary emission is dominated by the 
lines of [Fe II] (5.34 /im) in the IRAC 5.8 /um image, 
[Ar II] (6.99 /im) in the IRAC 8.0 /im image, [O IV] 
(25.89 /im), [Fe II] (25.99 /im) and dust emission in the 
MIPS 24 /xm residual image, and mostly dust emission 
in the MIPS 70 /im image with some contribution from 
[O I] (63.18 /im) and [O III] (88.36 (im) lines. We inte- 
grated the IRS spectra extracted from the synchrotron- 
subtracted spectral cube at the three positions in Fig- 




Fig. 4. — Map of the synchrotron-subtracted line-free region be- 
tween 27-32 (an, representing continuum emission from dust in the 
LL slits (magenta) overlaid on the MIPS 24 (im image (blue). The 
total coverage of the LL slits is shown in Figure [l] The south- 
ern and northern edges of the LL slits are marked by white dashed 
lines. The brightest dust emission is seen to spatially coincide with 
the brightest filaments. The extraction regions for the dust spectra 
in Figure [5] are overlaid in black and numbered. 

ure 0] over the MIPS 24 bandpass. We repeated 
this integration for line-subtracted spectra in order to 
compare the values and estimate the contribution of line 
emission in the MIPS 24 /im image. We find that approx- 
imately 27% of the residual emission comes from lines in 
Region 1, 54% in Region 2, and 48% in Region 3. This 
suggests that on average, about half of the emission in 
the MIPS 24 /im synchrotron-subtracted image comes 
from dust. 

The IRS spectral map of the dust continuum over- 
laid on the MIPS 24 /xm image is shown in Figure |4j 
The dust emission spatially correlates with the ejecta fil- 
aments, with the brightest dust emission concentrated 
along the brightest filaments. This is consistent with 
previous optical studies of the Crab Nebula that found a 
correl ation of dust absorption features with the filament 
cores (jFesen fe Blairlll99(H; [Blair et alJU997t iHester etaH 
ll99Ch . 

5. DUST COMPOSITION AND TEMPERATURE 

The extracted IRS spectra are shown in Figure [5] and 
they clearly show a rising dust continuum above the syn- 
chrotron power law. The residual dust emission shows 
no evidence for PAH features or for any obvious silicate 
emission features, other than those caused by the extinc- 
tion along the line of sight. The dust spectrum is fea- 
tureless and shows no emission features typical of ejecta 
dust in other young remnants (|Rho et al.ll2009D . 

In order to model the residual dust emission, we 
first corrected the spectrum for extinction along the 
line of sight us i ng th e mid-IR extinction curve from 
IChiar fc Tielensl (|2006T >. The extinction can be observed 
as the dip in the synchrotron power-law spectrum around 
10 /urn. We then fitted the synchrotron-subtracted spec- 
tra with a dust emission model u sing two different g rain 
composit ions; amorphous carbon (iZubko et al.ll2004T ) and 
silicates ([Weingartner fc Draind 120011) . Both composi- 




10' 1 — ■ — 11, 1 — 1 — 

10 30 90 

Fig. 5. — Model fits to the spectral continuum and the 70 fim 
data point for three extraction regions shown in Figure [4] Spectra 
from regions 1, 2, 3, and 4 are shown in plots (a), (b), (c), and 
(d), respectively. The model is silicate dust (dashed line), plus a 
synchrotron component whose index and normalization have been 
fixed to the values derived from the IRAC data (dotted line). The 
red line represents the best fit two-component model corresponding 
to a best-fit grain temperature of 53 K, 55 K, and 57 K, for plots 
(a), (b), and (c). There is an additional dust component with a 
temperature of ~ 136 K in Region 1. The amorphous carbon grains 
provide equally good fits to the data. The best-fit values for both 
compositions are listed in Table [3] For comparison, the spectrum 
in plot (d) is from a region centered on the PWN that does not 
show a significant excess from dust. The red line in plot (d) is the 
synchrotron component only, extrapolated from the 3.6 fim image. 

tions produced equally good fits to the observed spectra 
and the 70 /im data point. The silicate composition re- 
quires an average temperature of 55 ± 4 K, and amor- 
phous carbon grains an average temperature of 60 ± 7 
K. We do not see significant changes in temperature be- 
tween the three regions. 

There appears to be some excess emission between 20- 
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TABLE 3 
Dust Temperature and Mass 



Composition Best-Fit Temperature (K) Average Dust Total Dust Mass in 

Region 1 Region 2 Region 3 Temperature (K) the Crab Nebula (Mq) 

Silicates 53 ± 3 (136 ±12) 55 ± 2 57 ± 2 55 ±4 $-*ti%) x !0" 3 

Carbon 57 ± 6 (119 ±12) 61 ± 3 63 ± 2 60 ± 7 (3.2l^) x 1CT 3 

Note. — The temperature uncertainties in the three regions represent 1.6 a statistical uncertainties 
from the fits. The total mass values have been calculated by assuming that all of the dust in the Crab 
Nebula emitting at 24 /im is at the average observed dust temperature. In Region 1, the temperatures in 
parentheses are best-fit temperatures for the hot dust component that comprises less than 1% of the total 
mass in that region. 



30 jum in regions 1 and 2 that suggests that an additional 
warmer dust component with a much smaller mass may 
be contributing to the spectrum. While the addition of 
this second component did not significantly improve the 
fit for Region 2, an additional component with a tem- 
perature of 136 ± 12 K for silicates, and 119 ± 12 K for 
carbon grains did improve the fit in Region 1. If present, 
this warmer component represents only a tiny fraction of 
the total dust mass, on the order of 0.2 %. Such emis- 
sion may be explained by stochastically heated smaller 
grains. Since the spectrum below 15 /im is dominated 
by synchrotron emission, the SL data were not sensitive 
enough to constrain the emission from PAHs and the 
small, stochastically heated grains. 

The best fit dust temperatures for each region and each 
grain composition are listed in Table [3] The best fit dust 
models are shown in Figure [5] As described in Section 
13. H the power-law component that is added to the mod- 
els in Figure [5] was derived from the synchrotron index 
map in Figure [2] and the synchrotron-dominated IRAC 
3.6 fim image. The spectral index a (normalization) in 
the derived power-law component shown in Figure [5] is 
0.4 (15.0 MJy/sr) in Region 1, 0.6 (11.8 MJy/sr) in Re- 
gion 2, 0.7 (9.7 MJy/sr) in Region 3, and 0.36 (21.9 
MJy/sr) in Region 4. The normalization value is the 
surface brightness of the synchrotron component at 3.6 
/xm measured in each extraction region. 

6. DUST HEATING 

The dominant heating source for the dust in the Crab 
Nebula is synchrotron radiation from the pulsar wind, 
while collisiona l heating by the gas in t he filaments plays 
a smaller role (jDwek k, Werner! 1 19811 ). In order to de- 
termine if the PWN can heat the dust to the observed 
temperature, we calculated the expected dust heating 
rate by the PWN in the Crab Nebula and compared it 
to the cooling rate of dust at the observed temperatures. 
We assume that the filaments are optically thin to the in- 
coming radiation. A rough estimate of the optical depth 
gives t < 1, assuming a filament size of ~ 10 16 cm, a 
density of 2000 cm" 3 and a dust-to-gas mass ratio of 0.1 
(|Sankrit et alj fl998). Because of the low r, we neglected 
any internal absorption in the nebula. The radiation ab- 
sorbed by a dust grain is then given by the heating rate 

H = Aira 2 J L v Q{v, a)dv 

where L v is the non-thermal specific luminosity of the 
Crab Nebula's PWN at a given frequency, a is the grain 
size, Q(v, a) is the frequency and grain size dependent 
absorption coefficient for a given grain composition, and 



r is the distance between the radiation source, i.e. the 
pulsar, and the dust grain, assumed to be between 57"- 
100", or 0.55-1.0 pc for a distance of 2 kpc. This distance 
range was chosen based on the location of the ejecta 
filaments in 3-dim ensional models of the Crab Nebula 
(jCadez et all I2001 . The non-thermal br oadband spec - 
trum of the Crab Nebula summarized in iHesterl (2008) 
was used for the heating luminosity L v . In calculating 
the heating rate of the dust, we integrated the product of 
L v and the dust absorption efficiencies up to an energy 
of about 0.6 keV. The fraction of the energy deposited in 
the dust at higher energies makes a negligible contribu- 
tion of less than 1 % to the heating of the dust because 
of the combined effects of decreasing L v with energy, and 
the decreasing efficiency of the energy deposited by the 
photoelectrons in the dust (jDwek fc Smith 19961). 
The radiative cooling rate for a dust grain is given by 

L gr = Ana 2 J nB v (T)Q(u, a)dv, (4) 

where B v is t he Planck function and T is the dust grai n 
temperature (|Dwek fc Werner! Il981t iDwek et al.ll2008f ). 
In equilibrium, the heating rate H (r, a) and the cool- 
ing rate L gr (T,a) should balance each other. Since the 
distance from the heating source r and the dust temper- 
atures T are both measured parameters, and since both 
rates also depend on the grain radius a, we can calculate 
a grain size for which H{r, a) = L gr (T, a). This quantity 
then represents the grain size that is required in order 
for the PWN to heat the dust to the observed dust tem- 
peratures in Table [3] 

Due to the uncertainties in the fitted dust tempera- 
tures and the spread in the possible distance to the heat- 
ing source that ranges from 0.55-1.0 pc, there is actually 
a range of possible dust grain radii for each grain com- 
position. This is demonstrated in Figure |5] that shows 
the heating and cooling rates as a function of grain ra- 
dius for silicate and amorphous carbon grains. The solid 
black lines represent the heating rate at the two distance 
extremes, while the gray band represents the cooling rate 
for the spread of best-fit dust temperatures. The overlap- 
ping regions of the heating and cooling bands represent 
the possible dust grain sizes. 

It can be seen from Figure |6] that the silicate dust 
would need to be located at the closest distance extreme 
from the heating source in order to be heated to even 
the low-end of the observed temperatures. Since the fil- 
aments are distributed at some range of distances from 
the center of the Crab Nebula's PWN, the grain com- 
position may mostly be carbonaceous, or at least a mix- 
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Fig. 6. — The heating rate of dust by the PWN's broadband synchrotron spectrum, as a function of grain size, for a source distance 
ranging from 0.6-1.0 pc ( Cadcz ct al. 2004) is shown by the solid curves. The cooling rate of dust at the observed, best-fit temperature is 
shown by the gray band, where the width of the band represents the uncertainties in the temperature measurements listed in Table [3] The 
place where the heating and cooling rates cross gives the expected grain radius of the emitting dust. 



ture of silicate and carbonaceous material. Theoretical 
models based on classical nucleation and grain growth 
theory suggest that carbon dust is the most abundant 
dust species in Type IIP SNe with lower mass progen- 
itors (Ko zasa et al.l 120091 ) . so a carbonaceous composi- 
tion would not be unusual for the Crab Nebula, which 
is thought to have been produced in a Type IIP explo- 
sion with a 8-12 M w progenitor (e.g.. [Davidson fe Feseiil 
1985; IChevalieri[200alMacAlpine fc Satterfieldll2008j) ~ 

The most interesting thing to notice in Figure |5] is that 
the required grain radii are fairly small for both silicate 
and carbon grain compositions. The required grain radii 
are < 0.015 fim for silicate grains, and < 0.05 fim for 
amorphous carbon grains. Dust formation takes place 
in the rapidly expanding SN ejecta that are internally 
heated by radioactivties.The composition and sizes of the 
grains depend in a complex way on the kinematics, ther- 
modynamic history, and composition of the ejecta, as 
well as the initial number of nucleation centers around 
which the dust grows. The studies of dust properties in 
SNRs can therefore be used to relate the inferred grain 
sizes to the complex physical processes in the ejecta and 
provide some insight into the properties of the progenitor 
star. 

iKozasa et all (|2009D and lNozawa et al.l (|2010l ) have in- 
vestigated the dependence of dust formation in SN ejecta 
on the type of core collapse SN. They compared the mass 
and properties of dust formed in Type IIP SNe with 
massive H-envelopes and Type lib SNe with low mass 
H-envelopes and found that the total dust mass of 0.1- 
0.7 M Q does not depend on the thickness of this outer 
envelope, but that the grain size is strongly dependent 
on its thickness. A less massive envelope allows for a 
much higher expansion velocity of the He-core, leading 
to a much lower ejecta gas den sity, too low for the l arger 
grains to form. For example, IKozasa et al.l (2009) find 
that the radii of dust grains for a 0.08 M H-envelope are 
below 0.006 /im, while the radii of grains formed in a SN 
with a 10 Mq H-envelope are larger than 0.03 /im. The 
low initial mass of the Crab Nebula's progenitor and the 
relatively low expansion velocity of the optical filaments 
suggest that the Crab Nebula is a result of a Type IIP 
SN re.g. lChevalie"nl2005tlMacAlDine fc Satterfieldl 12008') . 



The models would therefore predict somewhat larger 
grain radii then inferred from our observations, unless 
the Crab's progenitor experienced substantial mass loss 
prior to the SN event. 

A low dust grain size also has important implications 
for the survival of SN dust and its injection into the 
ISM. The majority of dust grains with small radii are 
expecte d to be destr o yed by sputtering in the SN reverse 



shock dDwekl 120051 : IKozasa et all 120091: iBianchi et all 



120091 : iNozawa et al.l I2010T ) . Hydrodynamic simulations 
of grain destruction by SNR reverse shocks show that 
grains with initial radii below 0.1 / xm are almost com- 
pletely destroyed ()Silvia et al.l l2012f ) . This implies that 
the majority of the dust observed in the Crab Nebula 
by Spitzer should eventually be destroyed in the reverse 
shock. However, it is possible that clumping of gas and 
dust in the Crab's filaments may provide some shielding 
for the dust grains. 



7. DUST MASS ESTIMATE 

In order to estimate the total mass of dust in the Crab 
Nebula, we first estimated the contribution of dust emis- 
sion to the total integrated flux. After subtracting the 
synchrotron emission extrapolated from the IRAC 3.6 
/jm image (Figure |3|), the total residual contribution from 
dust and line emission at 24 /xm is 24.3 ± 2.4 Jy, or about 
40% of total flux. As described in Section @J we calcu- 
lated the fraction of line emission in three regions shown 
in Figure |4] We find that the dust continuum emission 
accounts for 73%, 46%, and 52% of the total residual 
emission in regions 1, 2, and 3, respectively. In order 
to estimate the dust fraction over the entire nebula, we 
assume an average value of 57% of dust emission, which 
gives 13.9 Jy for the total dust emission in the nebula at 
24 /im. We use this value, and the average of the best-fit 
temperatures for the three regions (Table [3|) to approx- 
imate the total mass of dust in the Crab Nebula. We 
used the equation 



M dus 



F u (v)d 2 



(5) 
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where F v (y) is the total IR flux at 24 /im, n v (y) is the 
grain opacity equal to 3 , d is the distance, B v (i>, T) 
is the Planck function evaluated at the grain tempera- 
ture, and p is the grain density. We used the opacity 
values f or silicate and amor phous carbon grains summa- 
rized in iZubko et al.l ([2004) . The dust mass values for 
both grain compositions are listed in Table [3] The total 
estimated mass for silicate grains is (2.4±f;|) x KT 3 M , 
and the total mass for amorphous carbon grains is 
(3.2^2^) x 10~ 3 M Q . We note that these estimates as- 
sume a single average dust temperature over the entire 
nebula, and that the actual composition of the dust is 
likely to be a mixture of the two composition. We es- 
timate the total dust mass in the Crab Nebula to be 
(1.2 — 12) x 10~ 3 M Q depending on the assumed compo- 
sition. 

The derived mass values are well below the theoret- 
ical pre dictions for Type II SNe. Classical n ucleation 
theory (jNozawa et al.ll2003t iKozasa et al.H2009Q and the 
theory of chemical kinemati c approach for the forma- 
tion of molecular precursors HCherchneff fc Dwekl [20101 
predict that each SN should form 0.1-0.7 M Q of dust. 
Type IIP SNe should form on the order of 0.1 M Q of 
dust, at least an order of magnitude more than what 
we observe in the Crab Nebula. Furthermore, the small 
grain radii inferred from the observations indicate that 
much of the observed dust will eventually be destroyed 
by the SN reverse shock. It is possible that dust grains 
with a colder temperature reside outside of the visible 
Crab Nebula's radius. These grains would emit mostly 
at far-IR and sub-mm wavelengths, and observations 
with the Herschel Space Observatory and the Atacama 
Large Millimeter/sub- millimeter Array (ALMA) should 
provide a constraint on the possible colder dust mass 
component in the Crab nebula. 

8. CONCLUSIONS 

In this paper, we have analyzed the Spitzer Space 
Telescope low-resolution spectroscopy of several regions 
across the Crab Nebula with the goal of isolating the 
emission from SN produced dust and determining its 
physical properties and spatial distribution. The dust 
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emission is spatially coincident with the brightest fila- 
ments in the ejecta. The dust spectrum appears feature- 
less, with no obvious PAH emission features, or other 
dust emission features typical of ejecta dust recently ob- 
served in other SNRs. We fitted the dust emission spec- 
tra with dust emission models of two different grain com- 
positions; silicates and amorphous carbon grains. Both 
grain compositions provide equally good fits to the ob- 
served emission. The best-fit model yields an average 
temperatures of 55 ± 4 K for silicates, and 60 ± 7 K for 
amorphous carbon grains. An emission excess between 
20-30 /xm suggest that a very small amount of warmer 
dust is present in some regions, possibly as a result of 
stochastic heating of smaller grains. 

The dust in the Crab Nebula is predominately heated 
by the surrounding synchrotron radiation field, rather 
than collisionally heated by the gas. Our heating mod- 
els show that the dust grain radii need to be relatively 
small in order for the dust to be heated to the observed 
temperature. The grain radii in the Crab Nebula would 
need to be < 0.015 [im for silicates and < 0.05 /im for 
carbon grains, unlike what is expected in Type IIP SNe. 
Assuming that the dust is composed predominantly of 
silicates or carbon, the total dust mass produced in the 
ejecta of the Crab is estimated to be (1.2 — 12) x 10~ 3 M Q , 
at least an order of magnitude less than was is expected 
from theoretical predictions. Due to the small grain size, 
much of the dust may eventually be destroyed by sput- 
tering in the SN reverse shock. 

Existing Herschel observations and future observations 
with ALMA will provide constraints on the mass of colder 
and larger dust grains that may be present in the Crab 
Nebula. The orders of magnitude better sensitivity and 
resolution of the James Webb Space Telescope (JWST) 
will allow us to probe the smaller dust grains and PAHs 
in the Crab Nebula's filaments, which will provide further 
insight into dust grain destruction and processing. 
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